Self-consistent Coulomb effects and charge distribution of quantum dot arrays 
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This paper considers the self-consistent Coulomb interaction within arrays of self-assembled InAs 
quantum dots (QDs) which are embedded in a pn structure. Strong emphasis is being put on 
the statistical occupation of the electronic QD states which has to be solved self-consistently with 
the actual three-dimensional potential distribution. A model which is based on a Green's function 
formalism including screening effects is used to calculate the interaction of QD carriers within an 
array of QDs, where screening due to the inhomogeneous bulk charge distribution is taken into 
acount. We apply our model to simulate capacitance-voltage (CV) characteristics of a pn structure 
with embedded QDs. Different size distributions of QDs and ensembles of spatially perodic and 
randomly distributed arrays of QDs are investigated. 

PACS numbers: 73.63.Kv 
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I. INTRODUCTION 

The strong localization of few electrons in three di- 
mensions is one of the key features of quantum dots Q 
(QDs). This results in a variety of properties making 
them potential candidates for new semiconductor appli- 
cations such as laser and memory devices 010. One 
of those properties is the singular density of states, which 
relies on the alignment of the energy levels in the whole 
QD ensemble. On the one hand this requires a fairly ho- 
mogeneous QD size distribution, on the other hand the 
statistical occupation of electronic QD levels changes 
the energy of levels of the same dot and all neighboring 
QDs due to Coulomb interaction. 

The amount of Coulomb repulsion between QD elec- 
trons, especially at different QDs, is expected to depend 
on the details of the complicated bulk charge distribu- 
tions of the QD device. Additionally, the microscopic 
charging and discharging processes of the QDs depend on 
the energy differences between the QD levels and the non- 
equilibrium quasi-Fermi level, resulting in a very compli- 
cated self-consistent problem. The microscopic charge 
distribution in the QD array is never stationary due to 
the statistical character of microscopic processes, but for 
large ensembles of QDs or in a long time average a sta- 
tionary mean QD sheet charge density is attained. Stan- 
dard approaches use quasi one-dimensional models where 
the charge in the QDs is approximated by a homoge- 
neous sheet charge ja, • This is suitable to obtain the 
QD energy levels by a fit to experimental capacitance- 
voltage (CV) data @, where a Gaussian distribution of 
QD energy levels was assumed. Nevertheless the inter- 
pretation of the fitted broadenings remains controversial, 
since they may result both from size fluctuations and in- 
homogeneous charging effects of the single QDs. 

The model used in this work calculates self-consistently 
the charge distribution in an array of QDs, with re- 
spect to the three-dimensional potential distribution un- 
der bias conditions. From an ensemble average of many 
QD charge configurations, the mean stationary QD sheet 
charge is calculated, to obtain macroscopic quantities 



like the depletion capacitance and the CV characteris- 
tics. The use of the electrostatic Green's function method 
makes this expensive self-consistent problem solvable and 
charging effects in the QD array can be investigated sep- 
arately from structural fluctuations. 



II. THEORY 

Since only a small but essential part of the structure 
is not translationally invariant in the lateral (x,y) direc- 
tions, but most of the device is inhomogeneous only in 
the growth direction (z) the problem of solving the three- 
dimensional non-linear Poisson's equation is divided up 
into calculating two one-dimensional potentials <3? and 
$ B Q and a suitable 3D Green's function. The poten- 
tial $ is introduced as the self-consistent solution of the 
non-linear one-dimensional Poisson's equation 

e a z [e(z)9 2 $(z)] = -p($, E Fa , E Fp ) - a QD X Qn(z) (1) 

and the current equations neglecting generation and re- 
combination in the active region 



d 2 j„[$,£ F „] =0 
d z j p [<S>,E Fp ] = 0, 



(2) 



where the QD charge is treated as a sheet charge 
(per unit area) averaged over the lateral directions. The 
bulk charge density p is a function of the local poten- 
tial <f>, and the electron and hole quasi-Fermi levels E Fn 
and E Fp . The electron current j n and the hole current 
j p are calculated within a drift-diffusion model. Ohmic 
contacts are used to obtain boundary conditions for $, 
E Fn and E Fp . For further details see Ref. @- Here, eo 
and e denote the absolute and relative permittivities of 
the semiconductor material, and xqd is the characteristic 
function 



XQd(z) = 1/h for z G QD layer 
XQdO) = else, 



(3) 



2 



where h is the height of the QDs, accounting for the 
localization in z direction. 

(pBQ j g p ^ en tial which is felt by a single QD charge, 
originating from the bulk charge if the QD charges would 
not exist. $ bc ^(zqd), where zqd is the position of the 
QDs, can be conceived as a background potential for the 
QDs. <& B Q is calculated as $ but without from eqn. 
©■ 

In this paper the intrinsic QD level energy E™^ of a 
state 11 — {n x ,n y ,n z ) in QD a (where n x , n y , n z are 
the respective quantum numbers) is measured from the 
conduction band edge. In order to specifically account 
for the size dependence of the level energies in InAs QDs 
we use a fit to the results of k -p calculations as shown 
in Fig. [2 

The charging energy E^ 1 originating from the 
Coulomb interaction of the QD electrons is treated as 
a first order perturbation. Thus the QD electron energy 
can be expressed as 



Eg? = E ca (z QB ) - e$ B Q(z QD ) + E^ + E* 



char 



(4) 



Here Eqo denotes the intrinsic conduction band edge, and 
e > is the elementary charge. The Coulomb charging 
energy E^ r can be expressed by the use of the capaci- 
tance matrix C 01 as 



e: 



char 



0u 



(5) 



Here pp v denotes the occupation numbers of the QD state 
v in QD (3 {jpp v = 1 if occupied, ppi, = if not occupied). 
^afi,(3u ensures that an electron does not interact with 
itself. The charging energies obtained by summing over 
[3 = a will be referred to as intradot charging energy, 
where the sum over (3 ^ a is referred to as interdot charg- 
ing energy in this paper. The elements of the inverse of 
the capacitance matrix (C -1 )^ can be expressed by the 

use of the Green's function G [llj : 

(CT 1 )^ J J |^M| 2 |^(^')| 2 G(r,r')d 3 r'd 3 r 
v v 

Here and ipf3 V symbolize the QD electron wavefunc- 
tions. Assuming that the QD electrons are distributed 
homogeneously within cylinder shaped QDs with height 
h and a diameter D a , neglecting the details of the QD 
wave functions leads to 



1 



V a Vp 



G(r,r')d 6 r'd 6 r, (7) 



v a v, } 



where V a and Vp denote the volumes of two QDs a and 
(3. In our model the charging energies are the same for all 
states li, v within a particular QD. A Green's function G 
including the static screening effects by free bulk charges 
can be found by defining a screening parameter / scr as 



/scr(z) 



Of) 

94> 



z)). 



(8) 



Then G is calculated as 0] 



A r + d z (lne(z)d z ) + 



/scr(z) 

e(z)e 



G(z,r,z') (9) 



1 



6((0,z')-(r,z)), 



where A r is the Laplace operator and r is the radial 
distance in the x — y plane. The solution of eqn.@ has a 
cylindrical symmetry, since / scr and e depend only on z. 
Note that G is the potential of a unit point-like charge 
positioned at z' , r' = 0. 

The QDs are occupied statistically according to the 
level energies £§? and the contact quasi-Fermi level Ep n 
with the occupation probability 



f 



1 + exp 



\ k B T J 



(10) 



where T is the temperature and ks denotes Boltzmann's 
constant. Note that every change in the occupation num- 
bers p QM changes all QD level energies. The stationary, 
mean QD sheet charge density cr^ D used in eqn. Q is 
found as a charge configuration average 



qd = z£j_y^ * 



(11) 



i— 1 a, fj, 



where i denotes Ni simulated realizations of configura- 
tions of occupied QD levels, and A is the quadratic con- 
tact cross section simulated. Periodic boundary condi- 
tions are used in the four parallel and four diagonal di- 
rections in the QD plane. Since \J~A 1 which takes values 
from 1 /jm to 3 pm in our simulations, is always larger 
than the maximum expected screening length for this 
structure, the influence of boundaries is negligible. In 
other words, a QD charge does not "see" itself in the 
periodically continued plane. The stationary depletion 
layer capacitance is determined as in Ref. by the elec- 
tric field in the depletion layer given by $. 



III. RESULTS 

The pn structure considered here consists of highly 
doped n + and p + GaAs contact regions, as well as a 
3 x 10 16 cm -3 n~-doped GaAs region. A layer of self- 
organized InAs QDs is positioned 400 nm from the pn 
interface. The structure is displayed in the inset of Fig. 
^ Similar pn structures have been investigated experi- 
mentally by CV spectroscopy and DLTS measurements 
[l3L Il4| . We simulate arrays of up to 30 by 30 QDs 
with a QD sheet density of Nqd — 10 10 cm~ 2 for various 
temperatures and ensembles of QDs. The average QD 
diameter is 12 nm. A QD height of 3 nm is used. 
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FIG. 1: Size dependence of the electron levels in In As QDs, 
according to a fit to k ■ p calculations The QD levels are 
labelled by the quantum numbers (n x , n y ,n z ). Inset: The pn 
structure under investigation with embedded self-organized 
quantum dots. 




FIG. 2: Calculated capacitance-voltage characteristics of the 
pn diode under investigation with periodically ordered and 
equally sized InAs quantum dots with a diameter of 12 nm 
and a height of 3 nm for various temperatures. Plots for 
300 K, 150 K, 35 K and 4 K are offset by -14, -7, +7 and 
+14 nF cm" 2 for clarity. 



peratures less than 150 K. The peak at U w-2.2 V orig- 
inates from the discharging of the excited states, which 
is split since piezoelectric fields break the degeneracy of 
the states (100) and (010) (see Fig. |TJ). At U w-3.6 V 
the two ground state electrons are discharged successively 
leading to a double peak. Here, the energy degeneracy is 
broken by the Coulomb charging effects. The Coulomb 
splitting has also been observed experimentally for differ- 
ent structures 0, 0, 0] . The peak structure is broad- 
ened with rising temperature, resulting in a plateau in 
the CV characteristics at 300 K. Our calculations use a 
charging energy of about 20 meV. From the ground state 
peak splitting in the 4 K curve of about 100 mV we can 
estimate a lever-arm factor of 200 meV/V which, when 
applied to the width of the ground state peaks, results 
in an energy broadening of 2.5 meV. This value coincides 
with the unscreened Coulomb potential (in units of en- 
ergy) of two electrons at a distance of 100 nm from a 
QD (2 x 1.25 meV). This means that the broadening at 
4 K originates from the charge fluctuations within the 
four nearest neighbor QDs, in particular if every second 
one is charged. At 300 K, where occupation is rather in- 
sensitive to the QD level spacing, the full plateau width 
of about 500 mV corresponds to an energy broadening 
of 100 meV. This value can be attributed to the broad- 
ening in InAs QDs of the Fermi distribution function of 
about 2 ksT at 300 K of the two ground state levels. 
Consequently, at large temperatures, the charging energy 
is given by the averaged, statistically distributed charge 
stored in the QD array. This is the reason why the CV 
characteristics shows a flat plateau at 300 K and is hardly 
distinguishable from that of a quantum well [l9l | . 

Additionally, the ground state peaks show a shift to 
lower reverse bias for rising temperatures, where the ex- 
cited states peaks do not show this feature. This shift 
can be explained as follows: When the ground states elec- 
trons are discharged, interdot charging energies become 
more important than for the excited states (explained 
later in detail). The interdot interaction depends on how 
the free electrons can screen the interdot Coulomb poten- 
tial. This ability rises with lower temperatures leading to 
shorter screening lengths. Therefore the charging ener- 
gies are lower at low temperatures, and a higher reverse 
bias has to be applied to discharge the ground state elec- 
trons. 



Ideal dot distribution 



B. Impact of structural fluctuations 



When a reverse bias U is applied to the diode, bulk 
carriers as well as QD states will be successively dis- 
charged ■§]. First we simulate the structure at 4 K, 35 K, 
70 K, 150 K and 300 K with periodically ordered QDs 
and with all QDs having the same size. The resulting 
CV characteristics are shown in Fig. [2J The main drop 
of the capacitance with increasing reverse bias is due to 
the increasing bulk depletion region ^J. On top of this 
background three peaks can clearly be observed for tem- 



Now we investigate the impact of structural fluctua- 
tions of the QD ensemble such as position disorder and 
size fluctuations. At first we distribute the equally sized 
QDs randomly instead of periodically in the lateral plane 
leaving the QD sheet density constant. Compared to the 
characteristics without fluctuations (solid lines in Fig. 0J) 
at 35 K and 70 K the ground state peaks for randomly 
distributed QDs are only slightly smeared out (dashed 
lines in Fig.[3Jl. Since the electrons in the excited states 
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FIG. 3: Calculated capacitance-voltage characteristics at 
35 K, 70 K and 300 K without structural quantum dot fluctua- 
tions (solid lines) , with randomly distributed but equally sized 
dots (dashed lines), with size fluctuations of ±17 % (dash- 
dotted line) and ±33 % (dotted lines). The curves for 35 K, 
70 K with size fluctuations and 300 K are offset by +7, -7 and 
-14 nF cm" 2 . 
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FIG. 4: Mean number of electrons per QD as a function of 
bias for 35 K, 70 K and 300 K without structural quantum 
dot fluctuations (solid lines), with randomly distributed but 
equally sized dots (dashed lines), with size fluctuations of 
±17 % (dash-dotted line) and ±33 % (dotted lines). The 
curves for 35 K, and 300 K are offset by +2 and -2 electrons 
per dot, respectively. 



can be thermally activated more easily due to fluctua- 
tions in the interdot Coulomb interaction, resulting in 
an inhomogeneous charging of these levels, the respec- 
tive peaks are more strongly broadened. 

Next we simulate an ensemble of periodically ordered 
QDs with random QD diameters between 10 and 14 nm 
at 70 K, corresponding to a fluctuation of ± 17 %. The re- 
sulting CV characteristic (see dash-dotted line in Fig. |3J 
shows a very broad ground state peak. Again, the ex- 
cited states are more sensitive to the fluctuations, i.e., 
the respective peak disappears almost completely. The 
increase of the fluctuations of the diameter to ±33 % (8 
- 16 nm) leads to a broad ground state plateau in the 
characteristics (dotted lines in Fig. |3J) ■ Such a plateau 
has been observed experimentally for similar pn struc- 
tures [TsLIT^ I. and a comparison with our results indicates 
that the size fluctuations in self-organized QDs is of this 
order of magnitude. This size fluctuation corresponds 
to an energy fluctuation of the ground state of about 
100 meV (Fig.^), which coincides with the typical value 
for the broadening of the ground state in the photolumi- 
nescence signal of InAs QDs 13] . For room temperature 
the CV characteristics with ±33 % size fluctuation shows 
almost no deviation from the characteristics without fluc- 
tuations. Here, thermal charge fluctuations dominate the 
behavior. 

As already mentioned, the QDs are discharged by ap- 
plying a reverse bias. In Fig. 0] the mean number of 
electrons per QD is displayed as a function of bias for 
35 K, 70 K and 300 K. At 35 K and 70 K the QDs are 
charged with three electrons at zero bias and are dis- 
charged monotonically in a step-like manner at biases 
where peaks are observed in the CV characteristics (solid 
lines). Structural fluctuations lead to a broadening of 



these curves (dashed, dash-dotted and dotted lines). Ad- 
ditionally, with size fluctuations of ±33 % some QDs can 
store up to four electrons, leading to a higher average 
QD occupation. For 300 K the QDs are charged with 
less electrons, because of the lower degree of degeneracy 
of bulk electrons in the contact region, and consequently 
the lower lying quasi- Fermi level at higher temperatures. 



C. Variation of the charging energy 

The mean charging energy per QD electron as shown in 
Fig-El s hows an interesting, non-monotonic behavior. At 
70 K without fluctuations (solid lines) the charging en- 
ergy is at first constant for U <-0.5 V and then rises until 
the excited state electrons are rapidly discharged. Be- 
tween U =-2.2 V and t/=-3.6 V the charging energy again 
rises, and then drops to zero when the residual ground 
state electrons are discharged. The continuous increase 
of the charging energies can be explained by consider- 
ing the bulk electrons in the vicinity of the QDs being 
depleted with reverse bias. This depletion increases the 
interdot charging energy, and therefore the total QD level 
energies. Fig. 03 shows that these negative and positive 
differential charging energies persist even for ±33 % size 
fluctuations at 70 K. Also for 300 K this behavior can be 
observed, but it vanishes when taking into account addi- 
tionally ±33 % size fluctuations. The mean height of the 
potential barrier c[<I>(zqd) — $ bq (zqd)] at the QD sheet 
as obtained from the one-dimensional Poisson's equation 
is also displayed in Fig. 03 It shows the big difference 
of the charging energies between one-dimensional and 
three-dimensional calculations, especially if the QDs are 
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FIG. 5: The mean charging energy per QD electron for 
70 K and 300 K without structural quantum dot fluctuations 
(solid lines) , with randomly distributed but equally sized dots 
(dashed lines), with size fluctuations of ±17 % (dash-dotted 
line) and ±33 % (dotted lines) as calculated from the 3D sim- 
ulations. The average potential barriers as calculated from 
the ID Poisson's equation are also shown (double dot-dashed 
line). The curves for 300 K are offset by ±50 meV. 



charged with more than two electrons. But nevertheless, 
it shows the same differential behavior, which has an in- 
hibitory effect on leakage and recombinations currents 
in QD devices, and therefore may provide an explana- 
tion for the observed bistabilities in similar pn structures 

To explain the origin of the negative differential charg- 
ing energies in more detail, the charging energies arc plot- 
ted in Fig.[5]as a function of bias, assuming that all QDs 
are charged with one electron at any bias and adding one 
extra electron into a single QD. We also differentiate be- 
tween the intradot charging energy (solid line) originat- 
ing from the second electron in a QD, and the interdot 
charging energy (dashed line) from the Coulomb repul- 
sion with all neighboring QDs. The ratio of the Green's 
function and the unscreened Green's function (/ scr = 0) 
is also displayed in the two insets of Fig.||)]for U=-1.5 V 
and ?7=-4.5 V as a function of the lateral position r and 
the vertical position z as contour plots. As shown in 
Section ff, the Green's function determines the capaci- 
tance matrix and therefore the charging energies. As the 
Green's function spreads out more and more, especially 
in the lateral directions, the charging energies due to in- 
terdot charging effects become larger. The effect on the 
intradot charging energy is rather small, since the size 
of the QDs is small compared to the distance of the free 
electron gas and the respective screening length. But by 
depletion of the free electrons with the applied reverse 
bias, the interdot charging rises. The slope of the in- 
terdot charging energy can clearly be identified with the 
slope of the mean charging energy of Fig. for biases 
where the QD charge is constant but non-zero. 
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FIG. 6: The intradot (solid line) and interdot charging energy 
(dashed line) as a function of bias at 70 K. Insets: The Green's 
functions (normalized to the unscreened Green's function) at 
U=-1.5 V (left) and U=-4.5 V (right) visualize the effect of the 
depletion of the free electron gas upon the charging energies. 



Fig. shows the QD charge distribution within the 
QD layer for U=-2.2 V, U=-2.8 V and U=-3.6 V and the 
energy distribution of the ground state (000) at U=0.0 V, 
U=-2.2 V, U=-2.8 V and [7=-3.6 V for 70 K and without 
structural fluctuations. The normalizations of the distri- 
butions have been chosen differently for clarity. The QD 
charge distribution is inhomogeneous when the QDs are 
discharged at U—-2.2 V and C/=-3.6 V, but it is homoge- 
neous at U=-2.8 V. Fig. shows that the initial ground 
state energy at U=0.0 V is lowered due to discharging 
one electron at U ~-2.2 V. Then the energy rises (up to 
U ~-2.8 V), due to the rising interdot charging energies 
as argued above. The discharging of the ground state 
electron then leads to a lowering of the QD electron en- 
ergies again. Overall, the QD energy shows a zig-zag 
behavior when sweeping the bias. 

The inhomogeneity of the QD charge leads also to a 
line broadening of the electron states. For U =-2.8 V and 
U=().Q V the charge distribution is homogeneous, and 
the broadening vanishes. As shown above the magnitude 
of the Coulomb interaction between the QDs depends on 
the bias (Fig.[f)J). Therefore the line broadening at U=- 
2.2 V is still small even for an inhomogeneous QD charge, 
since the interdot interaction is close to zero. The broad- 
ening becomes larger (about 2 meV) at J7=-3.6 V, since 
at this bias the charge distribution again is inhomoge- 
neous, and the interdot interaction is large (see Fig.[BJ). 



IV. CONCLUSION 

In conclusion, we have proposed a model to calculate 
the self-consistent three-dimensional charge distributions 
and charging energies within a QD device. We have 
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FIG. 7: The electron distribution within the quantum dot 
layer (top) and the lineshape of the ground state (000) for 
different reverse biases at 70 K (bottom). The discharging of 
the dots and the bias dependent Coulomb interaction lead to 
a nonmonotonic shift of the energy distribution. The charge 
inhomogeneity in the quantum dot layer leads to a broadening 
of the energy distribution at U=-3.6 V. 



shown that the charging energies depend on the elec- 
tronic structure, which is determined by doping concen- 
tration, bias and temperature. The effects of size fluc- 
tuations of the dots and the role of a spatially periodic 
or random distribution of the dots on the capacitance- 
voltage characteristics have been investigated. We have 
found that size fluctuations have a much stronger influ- 
ence on the broadening of the peak structure of the CV 
characteristics than the positional disorder of the QDs. 
Our detailed simulations have shown that the QD energy 
can be a strongly non-monotonic function of the applied 
bias. 
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